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Based on the stochastic gravity, we study the loop corrections to the scalar and tensor perturbations during 
inflation. Since the loop corrections to scalar perturbations suffer infrared (IR) divergence, we consider the IR 
regularization to obtain the finite value. We find that the loop corrections to the scalar perturbations are amplified 
by the e-folding; in other words there appear the logarithmic correction, just as discussed by M. Sloth et al. On 
the other hand, we find that the tensor perturbations do not suffer from infrared divergence. 



I. INTRODUCTION 

Inflation provides a natural framework explaining both the 
large-scale homogeneity of the universe and its small-scale ir- 
regularity. Despite its attractive aspects, there are still many 
unknowns about the inflation theory, since in most models in- 
flation takes place on an energy scale many orders of magni- 
tude higher than that accessible by accelerators. This is why it 
is necessary to learn all we can about this high energy regime 
from the signatures left by inflation in the present universe 

m. 

However, when we consider the power spectrum of the cur- 
vature perturbation C, only by linear analysis, many inflation 
models predict the same results, which are compatible with 
the observational data, although the fundamental theories are 
quiet different. To discriminate between different inflationary 
models, it is important to take into account nonlinear effects 
ll5l- [l6ll . In particular, the classical perturbation theory pre- 
dicts that when we consider most inflation models, the curva- 
ture perturbation which is directly related to the fluctuation 
of the t emp erature of CMB, is conserved in the superhorizon 
region lfl8l - l20ll . In this case, the primordial perturbation is es- 
sentially characterized by the behavior of the background in- 
flaton field near the time of horizon exit. Although this fact 
makes the computation of the generated primordial pertur- 
bations simple, it makes it difficult to discriminate different 
inflation models. That is why the non-local dependence on 
the evolution of the background scalar field has been studied 
among the nonlinear quantum effects such as the loop correc- 
tions |9, 10]. Despite their importance, it is difficult to com- 
pute these non-linear quantum effects. This is because they 
contain integrations regarding internal momenta i Tll - fTill . Fur- 
thermore, there are several types of nonlinear interactions that 
induce loop corrections, such as self-interaction of a scalar 
field and interaction between a matter field and a gravitational 
field. Depending on the interaction term, we find different 
loop correction behavior 

Stochastic gravity may be well-suited to computing loop 
corrections induced from interactions between a scalar field 



'Electronic address: 'yuko'at'gravity.phys. waseda.ac.jp | 
^Electronic address: maeda'afwaseda.jp 



and a gravitational field. Stochastic gravity was proposed as 
a means of discussing the behavior of the gravitational field 
on the sub-Planck scale, which is affected by quantum mat- 
ter fields |21-30]. From our naive expectation, on this energy 
scale, the quantum fluctuation of the matter field may domi- 
nate that of the gravitational field. Based on this insight, Mar- 
tin and Verdaguer have presented the evolution equation of 
the gravitational field, which is affected by a quantum scalar 
field L24J . The effect induced by the quantum matter field is 
evaluated by the so-called closed time path (CPT) formalism 
ifsH - llsll . We integrate the action over quantum scalar fields. 
As a result the evolution equation of the gravitational field is 
described by a Langevin-type equation, which is called the 
Einstein-Langevin equation. In general, we need great effort 
to compute the loop corrections. In stochastic gravity, how- 
ever, by focusing on non-linear interaction between a scalar 
field and the gravitational field which gives the leading con- 
tribution, we can compute such loop corrections much more 
easy. Hence, in this paper, using the Einstein-Langevin equa- 
tions formulated in 1 24], we evaluate one loop corrections in- 
duced by a quantum scalar field. 

In our previous work i36ll . we applied this formalism to the 
linear perturbations, especially to the curvature perturbation 
C, which is important to consider the imprint on observational 
data. We find that it reproduces the same results as the predic- 
tion obtained by the quantization of the gauge invariant vari- 
ables iIstI mi], except for the limited case. Only when the 
e-folding from the horizon crossing time to the end of infla- 
tion exceeds some critical value, does the Einstein-Langevin 
equation not g ive the same result as that of the gauge invariant 
variables [39]. Hence, we evaluate loop corrections to scalar 
and tensor perturbations assuming that the e-folding is smaller 
than the critical value. 

In general, loop corrections contain a divergent part. In a 
quantum field theory in Minkowski spacetime, the divergence 
usually appears on the high energy scale. To discuss finite 
and physical quantities, appropriate regularization and renor- 
malization are required. Apart from such an ultraviolet diver- 
gence, there may appear another divergence in de Sitter (or 
quasi de Sitter) spacetime. This infrared problem is important 
because if we introduce an infrared cut-off to obtain a finite 
value, which gives a logarithmic correction. Such a logarith- 
mic correction amplifies the perturbations. We also find that 
there is a crucial difference between the logarithmic correc- 
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tions in scalar and tensor perturbations, which is related to the 
infrared divergence. 

In this paper, we consider a minimally coupled single-field 



J 



inflation as a simple slow-roll inflation model, whose action is 
given by 



S[9, 



2k\ 



C/^abcd 
„ abcd'^ 



r 



PBR-^{9'''da^db(l) + 2Vm 



(1.1) 



where = SttGb is the bare gravitational constant. The 
subscript "B" represents the values of bare coupling constants. 
After we regularize divergent parts and renormalize them, we 
set the renormalized constants as a = l3 — A = Q for simplic- 
ity. We also represent the renormalized gravitational constant 
by = SttG. In order to characterize the slow-roll infla- 
tion, we adopt two slow-roll parameters: e = —H/H^ and 
r]v = V^^/k^V. As for the time variable, we use the confor- 
mal time, r, and represent the time derivative by a prime. 

The paper is organized as follows. In Sec. [Ill we briefly 
review the basic idea of stochastic gravity and consider the 
properties of the Einstein-Langevin equation, which describes 
the evolution of the gravitational field affected by quantum 
scalar fields. Then we consider perturbations of the Einstein- 
Langevin equation around an inflationary background space- 
time. In Sec. Uni we discuss how the loop corrections de- 
pend on the potential of the scalar field in diagrammatical 
language. This part is independent from the computation of 
the loop corrections in the later sections. In Sec. HV] we 
consider the perturbation of the Einstein-Langevin equation. 
Solving this perturbed equation, we can compute the primor- 
dial perturbations generated from the quantum fluctuation of 
the scalar field. The quantum fluctuation of the scalar field is 
represented by the stochastic variable. In Sec. [V] we compute 
the correlation function of the stochastic variables. Taking 
into account the results in Sec. |IV| and Sec. |V] we evaluate 
the loop corrections to the scalar and tensor perturbations. The 
conclusion and discussion foUow in Sec. IVIII 



II. STOCHASTIC GRAVITY 

First we shortly summarize the basic points of stochastic 
gravity and the Einstein-Langevin equation derived in [24J . It 
is a generalization of the semi-classical gravity theory. As- 
suming that quantum fluctuations of matter fields dominate 
that of the gravitational field, they quantize matter fields but 
treat gravity as a classical field. The fluctuations of the gravi- 
tational field induced through interaction with quantum matter 
fields are taken into account as stochastic variables. To discuss 
such gravitational field dynamics, the CTP formalism is use- 
ful. They derive the effective equation of motion based on the 
CTP functional technique applied to a system-environment in- 
teraction, more specifically, based on the influence functional 
formalism of Feynman and Vernon. It is worth while noting 
that this Langevin-type equation is well-suited not only to un- 



derstanding the properties of inflation and the origin of large- 
scale structures in the Universe but also to explaining the tran- 
sition from quantum fluctuations to classical seeds. In addi- 
tion to the ordinary Einstein-Hilbert action, this CTP effective 
action contains two specific terms, which describe the effects 
induced through interaction with quantum matter fields. One 
is a memory term, by which the equation of motion depends 
on the history of the gravitational field itself. The other is a 
stochastic source £^ab, which describes quantum fluctuation of 
a scalar field. The latter is obtained from the imaginary part 
of the effective action, and as such it cannot be interpreted 
as a conventional action. Indeed, there appear statistically 
weighted stochastic noises as a source for the gravitational 
field. Under the Gaussian approximation, this stochastic vari- 
able is characterized by the average value and the two-point 
correlation function: 

{^abixi)^c'd'{x2)) = Nabc'd'ixi, X2) , (2.1) 

where the bi-tensor Nabc'd'{xi, X2) is called a noise ker- 
nel, which represents quantum fluctuation of the energy- 
momentum tensor in a background spacetime, i.e., 

1 



Nabc'd' {xi,X2) = -Re.[Fabc'd' ixi,X2)] 
= Zil^abixi) - {fab{xi)),fabix2) 



{Tabix2)mg], 
(2.2) 



where {X, Y} = XY + YX, g is the metric of a background 
spacetime, and the bi-tensor Fabc'd' {x, y) is defined by 

Fabc'd'{xi,X2) = {fab(xi)fc'd'{x2))[g] 

-{fab{xi))[g] {t,a'ix2))[g].(2.3) 

The expectation value for the quantum scalar field is evaluated 
in the background spacetime g. 

Including the above-mentioned stochastic source of £^ab, the 
effective equation of motion for the gravitational field is writ- 
ten as 



G'"'[9 + 6g] 



ah 



(2.4) 



where 5g is the metric perturbation induced by quantum fluc- 
tuation of matter fields and stochastic source ^ab is character- 
ized by the average value and the two-point correlation func- 
tion Eq. dTT] ). 
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Note that this equation is the same as the semiclassical 
Einstein equation expect for a source term of stochastic vari- 
ables (_ab- Furthermore, the expectation value of the energy- 
momentum tensor includes a nonlocal effect as follows. It 
consists of three terms as 

-2 J d^y./^)H-''^%]{x,y)5gM + 0{5g^), 

(2.5) 

where the expectation value of T'(i)"'' and H°-^'^'^ are defined 
below (Eq. (12.61 ) and ( 12.8b ). The evolution equation for a 
scalar field depends on the gravitational field. As a result, 
the expectation value of energy-momentum tensor (Eq. ( I2.5l l) 
depends not only directly on the spacetime geometry but also 
indirectly through a scalar field. When we perturb a spacetime 
as {g + 5g), two different changes appear in the right hand side 
ofEq. ( 12.51 ). The second term in Eq. ( 12.51) represents the direct 
change, which is described by fluctuation of the gravitational 
field 5g as 

{f^'^'''m,5g]{x)) 



-{( 



2e 



2a2 



(2.6) 



where Stp'^ is defined in terms of the quantum fluctuation of 
the scalar field ip as follows: 



(2.7) 



We can neglect this term safely on the sub-Planck scale be- 
cause this term is smaller by the order of (kH)'^ than the pre- 
ceding term. To derive this expression, we have used the back- 
ground evolution equation for a scalar field. 

The third integral term in the r.h.s. of Eq. ( 12.51 ) represents 
the effect from the indirect change and is characterized by the 
dissipation kernel, which is given by 

Habc'd' {XI,X2) = ^^l6c'<i'(2;i'^2) +i?ltM'(2^1'2;2) (2.8) 



= ^lm[Fabc'd' {X1,X2)] 

where Sabc'd' {xi,X2) is defined by 



(2.9) 
(2.10) 



Sabc'd'{xi,X2) = {T*fat{xi)fc'd'{x2))[g]- (2.11) 

T* denotes that we take time ordering before we apply the 
derivative operators in the energy momentum tensor. As 
pointed out in ll24l] . only if the background spacetime g sat- 
isfies the semiclassical Einstein equation, is the gauge invari- 
ance of the Einstein-Langevin equation guaranteed. Hence, 
in this paper, to guarantee the gauge invariance, we assume 



the background spacetime satisfies the semiclassical Einstein 
equation. 

The Einstein-Langevin equation Eq. ( 12.41 ) contains two dif- 
ferent sources. One is a stochastic source ^ab, whose correla- 
tion function is given by the noise kernel. From the explicit 
form of a noise kernel Eq. ( 12.21 ). we find that £^ab represents 
the quantum fluctuation of the energy momentum tensor. The 
other is an expectation value of the energy momentum tensor 
in the perturbed spacetime {g + 5g), which includes a memory 
term. The integrand of a memory term consists of a dissipa- 
tion kernel and fluctuation of the gravitational field. To inves- 
tigate the evolution for fluctuation of the gravitational field, 
it is necessary to calculate the quantum correction of a scalar 
field and evaluate the noise kernel and the dissipation kernel. 
Note that the noise kernel and the dissipation kernel corre- 
spond to the contributions from internal lines or loops of the 
Feynman diagrams, which consist of propagators of a scalar 
field and do not include those of the gravitational field in our 
approach. 



III. GENERIC FEATURE FROM THE VERTEX 
OPERATORS 

Before we discuss loop corrections to the primordial per- 
turbations in detail, we consider the dependence of loop cor- 
rections on a potential of a scalar field diagrammatically. 

First we divide a scalar field (p into the classical part (f)ci and 
the part of small quantum fluctuation tp. Inflation is mainly 
driven by (pci- Expanding a potential of the scalar field V{(j)) 
around (pd, we can write 



oo ^ 
^ — ^ m I 



m=0 



where the coefficient a^") is defined by 



(3.1) 



(3.2) 



Taking into account that during inflation ipci changes on the 
Planck scale, we have normalized a*^™) by the Planck mass. 

Similarly, we also perturb the gravitational field in the total 
action. Expanding the total action Eq. dl.ll) in terms of the 
fluctuation of the gravitational field Sg and the fluctuation of 
the scalar field ip, we find that the following vertices appear: 



v^^'>Sgi;, a^^'>Sg'ilj^ , a^^^Sgtp^, a'^^'^SgiP^ 



.(3.3) 



Although the kinetic term also includes fluctuation terms such 
as Sgip and Sgif^, since in this section we are interested in in- 
formation about V{4)), which is obtained through the coeffi- 
cients a'-™', we do not pay attention to these terms. Of course, 
when we evaluate the loop corrections in the later sections, we 
take into account all the fluctuation terms. 

The vertex operators given by Eq. ( 13. 3t correspond to the 
vertex diagrams depicted in Fig. [T] The solid line represents 
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FIG. 1: Vertices 



the propagation of the scalar field ip. The equation of mo- 
tion for the field in the interaction picture is discussed in Sec. 
IV Al As we will see later, the propagator in the inflationary 
universe is proportional to H^. Hence, each solid line con- 
tributes as {kH)^. On the other hand, the wavy line represents 
the propagation of the gravitational field 6g. The coefficient 
a(") is a coupling constant of the interaction described by the 
vertex of Sg ip"^. Since the higher loop graphs are suppressed 
further by (kH)'^, we can discuss the quantum corrections by 
an iterative perturbation method. 

To consider the evolution equation of the gravitational field, 
we integrate out only the degree of freedom of a scalar field. 
This means that when we evaluate the effective action in the 
CTP formalism, the gravitational field is treated as a classi- 
cal external field. Hence the gravitational field contributes as 
the external line but not as the internal line, in the effective 
action. We represent the gravitational field by the wavy line. 
Taking into account this fact, we find the leading contribution 
to the effective action is given by the diagram depicted in Fig. 
|2](1). The amplitude of this leading diagram is proportional 

(1) (2) 

(kH) 




FIG. 2: Feynman diagrams. 



to {a^^^nHY; it depends both on the Hubble parameter H 
and on the first order derivative of the potential, a'^^. The 
contribution from this tree-level graph corresponds to quan- 
tum corrections in the linear perturbation analysis, and it has 
been given in our previous work ll36ll . In the next leading or- 
der depicted by Fig. |2](2), the amplitude is proportional to 
(a'^^^)^(Ki7)"'. It depends on the Hubble parameter H and the 
second order derivative of the potential, a^'^\ In this paper, 
we evaluate this leading contribution among loop corrections. 
We just give a few comments on further order contributions. 
As seen from Fig. |2] (3), there are two different diagrams. 



whose amplitudes are proportional to (kH)^. Figure |2] (3a) 
depends on the third-order derivative of the potential, a'-^K 
while Fig|2](3b) depends on the fourth-order derivative of the 
potential, a^"*'. Figure|2](3b) comes from the self-interaction 
of the scalar field. On the other hand, the other three graphs 
(Fig. |2](1), (2) and (3a)) are due to the interaction between the 
gravitational field and the scalar field (Sg4>™'{m = 1, 2, 3)). 
From our discussion here, the higher-order loop corrections, 
although they are suppressed by the Planck scale, make it pos- 
sible to know the more information about the potential of the 
scalar field. If we can detect these loop corrections, it helps 
to discriminate many inflation models, even if they cannot be 
distinguished only from the linear perturbation analysis. 



IV. PERTURBATION OF EINSTEIN-LANGEVIN 
EQUATION 

Next we discuss the behavior of the loop corrections, espe- 
cially those in the superhorizon region. We consider the time 
evolution of the leading loop corrections, which is depicted in 
Fig. |2](2), in the superhorizon region. To calculate the loop 
corrections to the scalar perturbations and the tensor perturba- 
tions, we adopt the following metric form: 

ds^ - -a\r){l + 2A^Yf^)dr^~2a^{T)^<^>j^Y^j^drdx^ 
W{t) (7,, + 2H^^le,j{k)Y^)dx'dx^ , (4.1) 

where a and 7^ are the scale factor and the metric of max- 
imally symmetric three space. The scalar perturbations are 
described by A and (k/K)^, which are the so-called lapse 
function and shift vector, respectively. H^^ is the tensor per- 
turbation. The scalar perturbations are expanded by a com- 
plete set of harmonic function Yf^{x), which satisfies 

{A + k^)Y}^{x) = . (4.2) 

Using these harmonic functions, we find the scalar compo- 
nents of vector variables are expanded by 

Y,k^-k-'yk\r (4-3) 

The tensor perturbations are expanded by the basis eij (k), 
which satisfies the transverse traceless condition, 7*^ eij (k) = 
and k'^Eij (k) = 0. Hence our variables are now Af^ , <i>^ and 

TK 

Since this gauge choice fixes both the time slicing and the 
spatial coordinate completely, all physical variables with this 
ansatz are gauge invariant. This choice of the time coordinate 
is called a flat slicing, because the spatial curvature vanishes 
in this slicing. 

Note that because of nonlinear perturbations, the tensor per- 
turbations are not decoupled from the scalar perturbations, 
and the quantum fluctuations of the scalar field may amplify 
not only the scalar perturbations but also the tensor perturba- 
tions. 
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A. Scalar perturbations 

First we consider loop corrections to the scalar perturba- 
tions. In particular, we discuss the evolution of the gauge- 
invariant variable (. We focus on proper nonUnear effects, 
and then we neglect the contributions from the product of the 
linear perturbations. Then ( is related to the density perturba- 
tion in a flat slicing {Sf = Sp/p) as 



(4.4) 



This variable C, is gauge-invariant and turns out to be a cur- 
vature perturbation in a uniform density slicing. In the clas- 
sical perturbation theory, the energy conservation law implies 
that this variable is conserved in a superhorizon region for a 
single-field inflation iflsl - Eoll . Q is directly related to a gravi- 
tational potential at the late stage of the universe and then to 
the observed CMB fluctuations. Hereafter, when we need not 
clarify the mode fc, we neglect the index of the momentum for 
the perturbed variables. 

The density perturbation in the present slicing is given by 

5T% ^ ~p5fY 

= 5goc{f'''{x))[g]+g,,{{f^'^'>^m,5g]{x)) 

-2 J d%^/^)H°'^''^[g]{x,y)Sg,4y) + 2^'} . 

(4.5) 

Since the background energy-momentum tensor is given by 



{T'''{x))[g]^g^''{T\ix))[g]^~pg 



Oa 



(4.6) 



the direct contribution from the gravitational field is described 
as 

(f (i)™[</.[5],<5g](x)) = + £ + 0{{KHf)}pAY . 

(4.7) 

With these two relations ( I4.5l l and ( I4.7l i. the density perturba- 
tion in flat slicing is written as 



3/ 



2e ^ 2.^ 
— TT-A + - (dp„ 
3 p 



(4.8) 



where we have defined the density perturbations of the 
stochastic source and of the memory term as follows: 



Spm = 



cl^xe 



ikx 



X iJooc'd' [g\{x, y)g''''' g'^'^' 5g^, f, (y) 
k X 



d-^xe 



-500 r(x) 



(4.9) 



The Hamiltonian constraint equation gives a relation be- 
tween the gauge invariant variable A and the density pertur- 
bation (5/ as 



zK'K) 2 



(4.10) 



Using it, we eliminate A in Eq. ( 14. St , and find 

(l - l)'^/ - + 5p„,) + 0((fc/J{)2) . (4.11) 

Hence, in superhorizon region, the two-point correlation func- 
tion for (5/ is expressed in terms of four correlation functions 
of 5p^ and i.e.. 



+ '^^Prnki'^)^P(.pi'^)) + i^P^k'^T)^Pr-^p('^)) 



+ {^P„rk^T)^Pmp{'r)) 



Here we have used the relation 

e 



V{t) 



B. Tensor perturbations 



(4.12) 



(4.13) 



In a similar way to the scalar perturbations, the transverse 
traceless part of the fluctuation of the energy-momentum ten- 
sor is given by 



STl 



TT 



(p^(^)+p7^J'))}e^(fc)yfc 



(4.14) 



where we have defined the transverse traceless part of the 
anisotropic pressure both for the memory term and the 
stochastic variable ^1 as 



ikx 



9'' / d'yV^) 



X Hkjc'd' [g] (x, y)g'''^'g'^'f'Sge'f' (y) 
p^fe'j{k) = J d'xe-'^-^ [9jke''ix)]TT 



TT 

(4.15) 



Taking into account that the transverse traceless part of Ein- 
stein tensor is written as 



SG\ 



1 



^^ = ^[dr + 2'Kdr + fc' ]HPe'j[k)Yj^ , (4.16) 



we find the transverse traceless part of the Einstein-Langevin 
equation as 

{d'^ + 2Jidr + e)H^l{r) 

^J,kir). (4.17) 

The l.h.s. of this equation is the same as the evolution equa- 
tion for linear perturbation. In contrast to the linear pertur- 
bation analysis, where the tensor perturbations are decoupled 
from the scalar perturbations, the non-linear interaction cou- 
ples these two modes. That is why in the rh.s. there appears 
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the influence of quantum fluctuations of a scalar field. A linear 
second-order differential equation with a source term is solved 
by the retarded Green function constructed from two indepen- 
dent general solutions. In the present case, since we are inter- 
ested in the tensor perturbations amplified by quantum scalar 
fields, we assume that the tensor perturbations were absent in 
the beginning of inflation. This gives {Ti) — as the ini- 
tial condition. The two independent general solutions for Eq. 
(14.171 1 are given by 



a{r) 



a{r) 



(4.18) 



(4.19) 



where v'^ = | + 3e and x — —kr. Hence, we find the solution 
forEq.( l4T7] i as 

dr'G,etk(T,r')J,fc(T'), 
where the retarded Green function is given by 



Grctk{T, t') 



Wk{,T') 



{r~r') 
(4.20) 



with 



dT 



(4.21) 



Substituting general solutions into these equations, we obtain 
the corresponding retarded Green function as 

Grct k(T, t') 



(T~T ). 

(4.22) 

Here we have used the formula for the Hankel functions: 

Hl')ix)^Hi'\x) ^ A. .(4.23) 

ax ax TTix 



Substituting these expressions into Eq. ( 14. 19b , we find the 
tensor perturbations amphfied by quantum scalar field as 



- / dr 



,a{T') 



a{T) 



r' Im[i7(i)(x)i/(2)(x')] J,fc(r') 



V. NOISE KERNEL 

The scalar perturbation C and the tensor perturbation iJ^*^ 
are given by the stochastic variable and the memory term. To 
evaluate the correlation functions for C, and it is nec- 

essary to compute quantum corrections for the scalar field, 
which are imprinted on the noise and dissipation kernels. It 
is expected that the contribution from the memory terms 5pm 
is smaller than that from the stochastic variable 6p(^ by the 
order of magnitude of the slow-roll parameters. The reason 
is as follows. The dissipation kernel is defined as two-point 
function of the energy-momentum tensor. As is seen from the 
definition of the contribution from the potential term is 
suppressed by the slow-roll parameters. Also, as summarized 
in Appendix B of our paper |36], the Green function scales as 
(_^)|slow-roll parameter! in the superhorizon region. Then, 
the time derivative of this Green function is suppressed by 
the slow-roll parameters. Taking into account that only the 
contribution in the superhorizon region can accumulate on the 
time integral of the memory term, we can see that the con- 
tribution from the memory term, which is proportional to the 
dissipation kernel, is suppressed by the slow-roll parameters. 
Neglecting the contribution from the memory term, we find 
the density perturbation 5 j only in terms of the density pertur- 
bation of the stochastic variable as 



f V 



(5.1) 



Similarly, the tensor perturbation i/^*' is given by the trans- 
verse traceless part of the anisotropic pressure of the stochas- 
tic variable as 



H 



(*)rr) - -— 

Tk^^> ^ fc2 



dx'[^ 



l+e 



(5.2) 

In this section, we shall evaluate the correlation functions of 
5p^ and pn^. In Appendix lAl we calculate these correlation 
functions from the noise kernel. They are expressed in terms 
of the Wightman Green function for the interaction picture 
field in momentum space, G'^{ti, T2 ). First we determine the 
Green function, and then we evaluate the correlation function 
of Sp^ and p-TT J. 



A. Propagator 



As mentioned before, to compute the correlation functions, 
it is necessary to determine the Wightman function in momen- 

^ lm[Hi'\x)Hi^Hx')]J^f^ix') SP^'^^' 

G+(ri, r2)^VM(ri)V;^,(r2), (5.3) 



(4.24) 



Here we have used the fact that the scale factor scales as 
a(T) cx |t|^(^+'^^ We have also omitted the contribution from 
the subhorizon region because the Hankel functions oscillate 
where x is larger than one. 



where ipf^k (t) is the mode function of a quantum scalar field, 
which satisfies the wave equation 

i^/,fc"(T) + 2Jfi^/,fc'(T) + + a^K^Vwj^fA^) = . 

(5.4) 
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We solve this equation under the slow-roll condition. Intro- 
ducing a new variable as ipk{T) = a(r) ijjf^kir), this equation 
is rewritten as 

V'fc M + [fc' - {2 - e - ?7v(3 - £)}M2]^fc(T) = , (5.5) 
where we have used the relation 



a^K^V = a^n'-pi 1 



2 2 ^ 



(5.6) 



On the sub-Planck scale, we can neglect the term whose 
magnitude is smaller by the order of (kH)^ than that of the 
leading term. We also ignore higher-order terms with respect 
to the slow-roll parameters. So we do not include the time 
evolution of slow-roll parameters. Under these assumptions, 
the equation for ip becomes 



■4>{x) 



1 - 



2 + 3{e-r]v) 



i'{x)=0, (5.7) 



where x = —kr, and we have used ~ — 1/[(1 — £)t]. The 
general solution is given by the Hankel functions as 



CH^p\x)+bHf\x) 



(5.8) 



where 0^ = 9/4 + 3(£ — ryy), with two arbitrary integration 
constants C and D. This implies 



■0fc(T) 



a(r) 



CH'^p\x)+DH\'\x 



/3 



(5.9) 



We assume that the mode functions should have the same form 
as in Minkowski spacetime, i.e.. 



1 



-ikTi 



(5.10) 



when the wavelength is much shorter than the horizon scale, 
i.e., at very early stage of the universe. This fact may be true 
in the present gauge rather than the comoving gauge. Then 
the mode function and the Wightman function in momentum 
space are given by 



■{2J3 + 1)_ 



a{T) 



Hf\x) 



(5.11) 



(5.12) 



Setting fli — 1, we give the scale factor a(T) by a(T) = 
(Ti/T) Using this fact, we rewrite the Wightman function 

as 



G+(ri, T2) = —[—I 

TT {xiX2)^ (TlT2 



L+e 



Hf\x,)Hl^\x2) 



(2) 



) {l-efH^H\p{x,)Hf\x2). 

(5.13) 



Here we have used the relation 

-2 /-I ^\2n/*2 



(5.14) 



To compute the correlation functions, it is sufficient to con- 
sider the evolution of the Wightman function in the superhori- 
zon region. The behavior of G^{ti, T"2)_in the superhorizon 
region is summarized in Appendix B in Il36ll . 



B. Scalar perturbations 

Once the Wightman function is determined, we can com- 
pute the correlation function of 5p^ from Eq. ( IA12b . 

(^p^fe(r)%p(r))(4) 

J ri,T2 = 

{(ai)-^(W^'-g- (fc-g)) +a['^V,^U^] 



8 

X 



X Re 



G+(ri, T2)G+_^|(ri, r^) 



(5.15) 



where the number of the superscript (4) represents the power 
of (kH). We put the momentum superscript on the par- 
tial derivative operator This means, for example, d^^ oper- 
ates only on the Wightman function with the momentum q, 
G+(ti, T2). It is convenient to divide this correlation func- 
tion into the subhorizon part Ish{T,k) and the superhorizon 
part /sp(r, fe), which are defined by 

/,b(T,fc) = / d^'q 

Jqe[J<, 00] 

X {(a2)-'(a,,^a,^'-'^ -q-{k-q))+ a^^^V^'n^} 
X Rc[g+(ti, T2)G+ (n, r2)l (5.16) 

L ^ ll^-qi J Ti,T2=T 

Jqelo, M] 

X {ia^)-'[d,^d^^^-'^~q-{k-q)) + a'^^'WfU^} 
X [{a^r^^d^p,,^-'^ -q-{k-q))+ ^'V^^'^^^j 

. (5.17) 



X Rc 



G+(ri, r2)G+_^|(ri, t^) 



First we discuss the subhorizon part, Ish{T,k). Since we con- 
sider only the superhorizon mode as the momentum of the 
external line (fc), k is much smaller than the horizon scale 
Di. This implies that the momentum of the internal line q 
in /sb(T, fc), which is larger than 5{, is much larger than the 
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external momentum, k. Hence, we can approximate |q — fc| as 
q. So Ish {T,k) depends only on t, and then 

/sb(r,fc) - /sb(T) - pfc'/sb(T) cx • (5.18) 

Here we have separated the scale invariant part of k~^. Since 
the remaining part, k^Ishir), must be a function of — fcr, we 
find that even the leading part of Ish{T,k) decays as {—kr)^. 
Hence we can neglect the contribution from the subhorizon 
region. 

Consequently, the correlation function of 5p^j^{t) is evalu- 
ated only in the superhorizon region as 



Substituting this result into Eq. ( 15.191 ). we obtain the correla- 
tion function of 5p^]f (t) as 



vr {hH{t)Y 



(4) 



(t^vV^^^^ x^^'^^-'^ S{k+p) 



8 fc3 

1 1 - (i/i/fc)-^!^-^) 1 / fc \3+4(£-w)- 

3 2(e - r/y) 3 viK; 

(5.23) 

If e > ryy, it diverges when we remove the cut-off Hi. 



{5p^k{T)6p^p{T)) 



(4) _ t 



d{k+q)hp{T,k). (5.19) 



As seen in Appendix B in 1361 . in the superhorizon region, the 
Wightman function (ti , T2) is approximated as 



Gfe (ri, T2) ~ {X1X2) 



2 0x02 

1 (Tir2)^-^ 

2 aia2 



. (5.20) 



Substituting this expression into the definition of I^p (t, k) and 
neglecting the sub-leading terms w.r.t. the slow-roll parame- 
ters, we obtain /sp(t, k) as 



/sp(T,fc) = {tIvV^'k'] - 



2 1 |r|2-4/3 



4 a(T)'' 



d3q 



?G[0, M] g3+2(e-r,v) |fc - q|3+2(=-w) ' 

(5.21) 



Here we encounter the so-called infrared divergence problem. 
In the long wave limit {q — > 0), the integrand is approximately 
^-[3+2(£-J7v)] xhen this integral could be divergent depend- 
ing on the signature of (e — 7]y)|,68i,l. This is the infrared (IR) 
problem, which sometimes appears in the quantum field the- 
ory in an inflationary universe. When we use the scale in- 
variant power spectrum, in general we find this divergence on 
the loop corrections. Here, introducing the cut off by the ini- 
tial horizon scale, we just neglect the effects from the long 
wave modes whose comoving lengths are larger than the ini- 
tial horizon scale aiHi. We tentatively discuss this IR problem 
in Sec lVIII and elaborate this problem in |l69'l. 

After introducing the cut off Hi and integrating over the 
internal momentum q, we obtain a finite result. Using the 
loop integral Eq. ( IB4I ). whose detailed derivation is given in 
Appendix C, we find /sp(T, fc) as 

hp{T,k) 

1 1 - (i?,/fc)-2('^-'?v) 1 , k \3+4(e-w) 



C. Tensor perturbations 

Next we calculate the correlation function of the transverse 
traceless part of the anisotropic pressure of the stochastic vari- 
able, pTTcf. (t), which is given by Eq.( |A17| i in Appendix B, as 



(P^ik (^1 ) (fc )Pnpi^2) e'i ip)) 



(4) 



4(aia2) 



S{k+p) I £q (q'- 



{kq 



|2x2 



xRe 



G+(ri, r2)G|^_g|(ri, t^) 



(5.24) 



For the tensor perturbations, we also divide the correlation 
functions into the subhorizon part Jsh{Ti;T2,k) and the su- 
perhorizon part Jsp{Ti,T2,k), which are defined as 



Jsb(Tl,T2,fc) 

EE - T2) 



(fq [q 



I qe[K2,oo\ 

xRe 



(fc-g) 

fc2 



2\ 2 



(t2 - Tl) 



d^q q^ - 



[k-q) 



2x2 



xRe 



(5.25) 



Jsp{Tl,T2,k) 
= e{Tl - T2) 



{kq 



(2x2 



xRe 



G+(ri, r2)G|^_^|(Ti, T2) 



0{t2-ti] 



qe[HiMi 



d'q iq'- 



(fc-g) 

fc2 



2.2 



xRe 



G+(ri, r2)G|^_g|(Ti, r2) 



(5.26) 



2(^s — rjv) 3\D{J J Note that to compute the correlation function for the ten- 

(5.22) sor perturbation, H^^j^ (t) at conformal time r, it is necessary 
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to consider the correlation function of jJTr^j^ for two differ- 
ent times Ti and T2. This is because, as seen from Eq. ( I5.2l i. 
the expression of H^^ includes the time integral. Therefore, 
there are two different comoving horizon scales correspond- 
ing to the different times ti and T2. For the same reason as 
in the case of the scalar perturbations, we have introduced the 
IR cut-off Hi. Nevertheless, we can see later, for the tensor 
perturbations, even if we remove the IR cut-off, the loop cor- 
rections remains finite. 

By virtue of the same argument as that presented in the 
scalar perturbations, Jsb(Ti,T2,fc) contains only the decay- 
ing modes. To show this, note that if either — fcri or —kT2 is 
larger than unity, it does not produce cumulative contributions 
because of the oscillation of the Hankel function in subhorizon 
region, as mentioned in Eq. ( |4.24| l. Hence, it is sufficient to 
consider only the case where both — kri and — fcr2 are smaller 
than unity. If ri > T2, then the inner momentum q is larger 
than IK2 — ^l/''2- Hence, this implies that q is larger than 
k. Approximating |fe — q| as q, we find that Jsb(Ti, T2,fc) con- 
tains only the decaying mode as we have shown in the scalar 
perturbations. The same discussion is valid also in case of 
T2 > Ti . Hence, in order to compute the correlation function 
of pTT^*', it is sufficient to consider only the contribution from 
the superhorizon region, Jsp (ti ,T2,k). 

As with the case of scalar perturbations, substituting the 
approximation of the Wightman function in the superhorizon 
region into the definition of Jsp(ri, T2,fc), we find the contri- 
bution from the superhorizon region as 



-'"sp(Tl,T2,fc) = 



1 inr^r-'^ 

4 (0102)2 



(T1-T2 



d\ (q^ 
1 



{k-qf 



q3+2{e-rtv) \q - k\^+2{e-7jv) 



9iT2-n) £qU 

Jqe[H.Mi] ^ 

1 



k^ 



q3+2{e-r,v) \q - k\^+2{e~nv) 



(5.27) 



This loop integral is given by Eq. 
implies 



in Appendix C. It 



Jspin, T2,k) 

Stt 
15 



^-^Htk 



{X2 - Xi) X^ 
-X2)x-'+'^-'^ 



2»7v„-l+4e-2r,v 
^2 



-4(^1^2)'"^ 



(5.28) 



where we have used H{t)'^x~''^^ = H^. To derive this re- 
lation, we have taken the limit of Hi 0. No divergence 
appears. It means that the cut-off for the infrared region is 
not necessary. It is interesting to note that scalar perturba- 
tions have the infrared divergence, but tensor perturbations do 



not suffer from the infrared problem. In Sec lVIII we discuss 
the origin of this infrared divergence and the reason why only 
tensor perturbations does not contain it. 

As a result, we obtain the correlation function for the tensor 
part of the anisotropic pressure as 



- (0102) ^ 5{k 



27r 



15 (0102)' 



S{k +p 



P) Jsp(Tl,T2,fc 
0{X2 



+d{xi - X2)X^ 



:ixiX2) 



2-nv 



(5.29) 

In our computation, we have neglected the loop corrections 
from the tensor perturbations, because the amplitude of power 
spectrum for the tensor perturbations is smaller than that for 
the scalar perturbations by order of the slow-roll parameter e. 



VI. LOOP CORRECTIONS TO THE CORRELATION 
FUNCTIONS 



As shown in Eqs. ( 15.11 ) and ( 15. 2t , the leading parts of the 
correlation function of the density perturbation in a flat-slicing 
5 f and of the tensor perturbation H^^ are determined by the 

stochastic variables (5p^ and pTr^'''. We then have calculated 
the correlation functions of Sp^ and pir^ for the noise ker- 
nel. Combining these results, we shall evaluate the correlation 



(t) 



functions of 5/ and H}p 



A. Scalar perturbations 

To evaluate the correlation function of the density perturba- 
tion in flat-slicing S / , focusing on the proper nonlinear effects, 
we neglect the contribution from the product of linear pertur- 
bations. This density perturbation is related to the curvature 
perturbation in uniform slicing by Eq. ( |4.4| i. Hence, once 
we find the loop corrections to Sf, we also obtain the loop 
corrections to (. The curvature perturbation C is proportional 
to the gravitational potential in the late time of the universe 
and it is directly related to the fluctuation of the temperature 
of CMB. That is why it is important for us to consider this 
gauge-invariant variable among scalar perturbations. 

From Eq. ( 15.11 ) and (15.231 ). the loop corrections to the cor- 
relation function of the density perturbation are given by 



2 fc3 
rl 



1 



■^krf'^^ 5{k + 
{k/H,f^^-'^'' 



P) 



-3 2(e-r,y) . 

where we have used the relation of H{t)'^x~'^^ 



(6.1) 



Hi. In 



our previous work |36], we showed that when we solve ex- 
plicitly the Einstein- Langevin equation [24], which includes 
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an iterative aspect, it disturbs the constant evolution of ( in 
superhorizon region. Since this is a problem of the way to 
quantize the scalar field and the gravitational field, the same 
effects may exist in the present loop corrections. Hence, we 
restrict our discussions to the case when rjv log x is smaller 
than unity, i.e., we assume that (— fcr)'*''^ ss 1. Then, we find 
the correlation function of the density perturbation S / as 



the part of the integrand is approximated as 

2 

3^ 



){5 



1 - (k/H,) 



2(£ - T^v) 



}■ 

(6.2) 



Taking into account Eq. ( 14. 4t . we obtain the loop corrections 
to the correlation function of the curvature perturbation in uni- 
form density slicing C as 



(Cfe(r)Cp(r))(^) 



8 fc3 



1 1- (fc/ff,)^''^'"'^ 



2(£- W) 



The final result depends on the initial Hubble horizon scale, 
"Ki = Hi, which is introduced to remove the infrared diver- 
gence. The case with 2\e — Tyy | log(fc/J{,:) < 1 is particularly 
interesting. This, in other words, corresponds to the case of 
Nk < l/2\e - I, where Nk ~ log(A:/J£i) is the e-folding 
from the beginning of inflation to the horizon crossing time. 
In this case, this correlation function is approximated as 



(Cfc(r)Cp(r))(^) 

TT {KHkY 



8 k^ 



Nk . (6.4) 



Note that there appears the logarithmic corrections. These re- 
sults imply that although the one-loop correction is suppressed 
by [nHk)'^ and is smaller by the order of the {KHkY than 
tree-level effects, it is amplified by the e-folding A^^ from the 
initial time to the horizon crossing time, which can become 
large contrary to the e-folding from the horizon crossing time 
to the end of the inflation. However note that this amplifica- 
tion is derived by introducing the IR cut-off and the obtained 
loop corrections significantly depend on the choice of the IR 
cut-off. 



B. Tensor perturbations 

The tensor perturbation H^'^ is related to the source term 

pTrf\ and the correlation function of pTr^*'' is given by Eq. 
( 15.291 1. Then, integrating over x — — fcr, we obtain the cor- 
relation function and the amplitude of the tensor perturbation, 
which could be amplified by the quantum effect of a scalar 
field. To integrate over x, it is helpful to note the asymp- 
totic behaviour of the Hankel function when the argument x 
is smaller than unity. As summarized in Appendix B in ll36ll . 



Im 



(6.5) 



Using this approximation, we find the loop corrections to the 
correlation function of the tensor perturbations as 



(i^^*^(r)e^(fc)i^^'^(r)e^») 



(*) 



\(4) 



135 P 



d{k+p) f - - 15 a;2+2''^' 
V 6 



(6.6) 



It is interesting to note that there is no dependence on the in- 
frared cut off Hi in the tensor perturbations. Furthermore, the 
tensor perturbations are divergence free in the superhorizon 
region. We shall discuss the reason in the next section. In 
our computation, we have neglected the loop corrections from 
the tensor perturbations, since the amplitude of the tensor per- 
turbations are smaller than that of the scalar perturbations by 
order of the slow-roll parameter e. 



(6.3) 



VII. DISCUSSIONS 

Using the Einstein-Langevin equation proposed in JH], we 
calculate the loop corrections to the scalar perturbations and 
the tensor perturbations, which are amplified through the non- 
linear interaction between the scalar field and the gravitational 
field. Here we discuss the origin of the amplification. 

When we consider the loop corrections in inflationary uni- 
verse, there are two different divergences. One is the ultravi- 
olet (UV) divergence. Since this divergence is originated by 
short wave modes, such divergence also appears in the quan- 
tum field theory in a Minkowski background. In inflation- 
ary spacetime, there exists another divergence, which is not 
found in Minkowski spacetime. This is the IR divergence. To 
avoid this IR divergence, we have introduced the cut-off at the 
initial Hubble horizon size. Then the amplitude of the one- 
loop corrections to the curvature perturbation ( is amplified 
by the e-folding from the initial time of inflation to the hori- 
zon crossing time, i.e., the logarithmic correction. If this is 
true, this amplification may make it possible to detect these 
loop corrections. Then it will be a great help to clarify the 
fundamental properties of an inflation model. 

So far we have several discussions about this logarithmic 
corrections due to IR divergence. Early works about this 
problem are done by Boyanovsky, de Vega and Sanchez ||64| - 
16711 . They calculated one loop corrections by light scalar 
and fermion fields to the inflaton potential, and also evalu- 
ated those by the gauge invariant curvature and tensor pertur- 
bations. They found that there appear the IR enhancements 
both in the scalar field corrections and curvature perturba- 
tions, while both fermion corrections and tensor perturbations 
do not exhibit IR divergences. Weinberg also pointed that 
the loop corrections to the primordial perturbations behave 
at most logarithmic Afterward Sloth considered the 

loop corrections to the fluctuation of the scalar field in flat- 
slicing ifTll [T2I1 . To avoid IR divergence, he introduced the 
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cut off by the initial horizon scale. As a result, he found that 
the loop correction is amplified by the e-folding from the ini- 
tial time of inflation to the horizon crossing time, which is 
also found in this paper from the analysis based on stochas- 
tic gravity. Following Sloth, Seery readdressed this problem 
more carefully lfl3i[l4ll . In particular, he analysed the evolu- 
tion in the superhorizon region using the 6N formula jTsi lssll . 
and improved his results. In this paper, we have computed 
the loop corrections by stochastic gravity, and found the simi- 
lar logarithmic corrections for scalar perturbations. The same 
logarithmic behaviours have been found in other interacting 
systems [59-61]. However, we should note that the IR prob- 
lem requires the careful treatment and the way to evaluate IR 
effects is controversial (6^j6^. 

Recently, Lyth has claimed that, to avoid the assumptions 
on unknown parts of the universe, the calculation about loop 
corrections should be done inside a comoving box, whose size 
L is not too much bigger than the present horizon scale ifTsll . 
The IR corrections are significantly reduced, although we still 
find the logarithmic behaviours. Furthermore, Bartolo et al. 
claimed that a stochastic approach plays a crucial role to deal 
with this problem IIITH . In relation to their insists, we should 
stress our interesting results. That is, although the scalar per- 
turbations are amplified by the logarithmic corrections, the 
tensor perturbations are not. Even if we remove the IR cut- 
off, the IR divergence does not appear in tensor perturbations. 
This difference between the scalar and tensor perturbations 
seems to be related to the origin of these logarithmic correc- 
tions. 

To consider the origin of this logarithmic corrections due 
to the IR cut-off, we first mention the prediction in stochas- 
tic inflation ll40l - l42[ I4i - l5lll . In stochastic inflation, the long 
wave mode k with k < aH of the scalar field couples to the 
short wave mode k with k > aH through the nonlinear self- 
interaction of the scalar field. Then the long wave modes are 
affected by the quantum fluctuation of the short wave modes. 
As a result, the long wave modes come to show stochastic 
behavior This stochastic behavior of the long wave modes 
affects the background quantities. 

In our case, due to the nonlinear interaction between the 
gravitational field and the scalar field, the long wave modes 
and the background quantities come to show the stochastic 
behavior Since the scalar perturbations are defined as the de- 
viation from the background quantities, the stochastic fluc- 
tuation of the background quantities affect the behavior of the 
perturbed variables. As a result, it induces the logarithmic sec- 
ular evolution of the perturbed variables. On the other hand. 



there are no background tensor modes, and then the one loop 
corrections to the tensor perturbation can avoid being affected 
by the background stochastic fluctuations. Furthermore, al- 
though in this paper, as the simplest step for treating the IR 
divergence, we have simply neglected the long wave modes 
with —k'Ki > 1, the infrared modes require the more care- 
ful treatment. We will propose the way to regularize the IR 
corrections in ll69ll . 

There is another notable difference between the scalar and 
the tensor perturbations. As pointed out in our previous work, 
in the present approach of stochastic gravity the longitudinal 
part of gravitational field is included iteratively. This affects 
the behavior of perturbations in the superhorizon region. In 
particular, the curvature perturbation deviates from constant 
when the e-folding from the horizon crossing time exceeds the 
definite critical value ( |slow-roll parameter] ^^). Since this is 
the problem of the way to quantize the gravitational field and 
the matter field, the loop corrections to the scalar perturbations 
are also influenced by the nonexistence of the longitudinal part 
of the quantized gravitational field. In fact, as shown in Eq. 
( I6.lt . the one loop correction to the curvature perturbation C 
evolves as x*^^' in the superhorizon region. On the other hand, 
as shown in Eq. ( 16.6b . the tensor perturbations do not decay. 
This means that even if we use the Einstein-Langevin equation 
in the present iterative way, it does not affect the one loop 
correction to the tensor perturbations. 
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Appendix A: Computations of tlie Noise kernel 



In AppendixlAl we calculate the noise kernel, which is defined by Eq. ( I2.2l i and (I2.3l l. As shown in SeclVIl we have to compute 
the correlation function of the density perturbation and the transverse traceless part of the anisotropic pressure of the stochastic 
variable ^ab- Then, we compute the (0, 0, 0', 0') component and the transverse traceless part of the (i, j, k', I') component of the 
noise kernel, i.e., -Fq^o' (^^i, 2:2) and F^j'^^, X2). Note that the noise kernel is computed from the quantum fluctuation of the 
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scalar field on the background spacetime. Decomposing the scalar field into (j) ^ ip, the energy momentum tensor 



Tab = V,0Vb0 - i5ab[V,0V^0 + 2Vi4')] 



(Al) 



is expressed as the classical part and the fluctuation part as follows: 

Tab = ri"" + STab ; 



ab 

Jl'-^ EE 6"j"^a'^<pli + ^gab^'ll - 9abV{<j)cl) 



T 



(A2) 



1 °° (m) 

STab = {SlVbi' + S'iVay^ + llabV^i>)a^cl + ^ a^N bi> - -gab^c^jV^^J - gabV{(bcl) V ^(kV-)"' • (A3) 

^ — ' to! 



The noise kernel, which represents the fluctuation of the energy-momentum tensor, can be expressed in terms of STab- In 
fact, substituting this decomposed energy-momentum tensor into the definition of Fabc'd' {xi , X2), we can express the two-point 
function Fabdd' , 2:2) as 

Fabc'd' {X1,X2) 

= {5T ab{xi)5T c' d' {X2)) - {5Tab{xi)){5Tc'd'{x2)) 

+ Tic'd'{5\G% + 5\G%,) + Tiab{5''dG%^. + CG%d) + r^abric'd'G%^,] 

-t- + U.^^/Lr.bg, — (^aij rjab^-ec'^-d' ~>>"2j Vc'd'^:af'^-b "I ^ VabVC d' ^-ef ^ 



(a2)^»?c'd'V'c(,; 



(aia2)^77Qb77c'd'V'ci,iVd,: 



~(2)~(2) 

aW4''.^G^ + ^^ 
2! 



(A4) 



where G^ = {^l,\^l;H{xl)1pH{x2)\^) is the Wightman Green function for the interacting system, which is defined as the two- 
point function of the Heisenberg field ^Ijh- Here we have redefined the coefficients ci™, including the divergent part G^(x, x) 
as follows: 



(3) (5) 2 



5(2) ^ ^(2) 



J4) 



(A5) 



Strictly speaking, we have to renormalize these divergent terms into the coefficients of the potential V{(j>). In this paper, we 
assume that these divergent parts are removed by an appropriate renormalization procedure. So the finite part of these radiative 
correction terms is much smaller than the leading term among the coefficients a^"^\ Here we neglect them and approximate 

as a^™). 

In our previous work ll36ll . we discussed the linear perturbations which are proportional to {kH)^. In this paper, we consider 
the leading loop corrections which are proportional to {kH}'^. The self-interaction part contributes to the effective action from 
the order (kH)^, which is depicted in Fig. |2](3). In these diagrams, the solid line represents the interacting picture field which 
satisfies the equation 



[9^0 + iD- 2)J{do - V" + a'V"{<f>,i)] = 



(A6) 



When we compute the effective action up to the order of {kHY, we can replace the Wightman function for the Heisenberg field 
G^ {xi,X2) to the Wightman function for the interaction picture field G+(xi,2:2). Then, the parts of Fabc'd' {xi,X2), whose 
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orders are (kH)^ and (kH)*, are given by 

- {aifvabVcisCi2ku2 G+ + 5%, G+ + r;,,^, G+ ) 

- (a2)'?7c'd'^^"c/.2ai0c/,ia2'^ G+ + ^Ofc G+ + r;„b G+) 

+ (aia2)^77a6'7c'd'Vc/,iK/,2 o'l^a^^K^G'^ (A7) 

and 

J^i^cM' (^1,2^2) = G+,,G+,, + G+,,G+,, - (ai)2r;,,G+,,G^%, - (a2)'ry,,,,G+^,G+ + I^i|^,y,,^,,,,G+^,G+^^/' 

- {aifi^abV.i^ia^l^n' { G+ G+ - i (a2)2rye,d'G+ G+-/' } 

- (a2)'r7c'd'Ki,2a^'^K' { G+G+ - ^ (ai)2,?,bG+G+-^ } 

(2) (2) 

+ {aia2fiiabVc'd'VcLiVcL2 («;2G+)2 , (A8) 

respectively, where the superscript numbers represent the powers of (kH). To compute the correlation functions of the primordial 
perturbations in momentum space, it is convenient to use the Fourier-transformed Wightman function given by 

G+{r,,T2)^yj^f,in)rfj,{r2), (A9) 

where the mode function tpj: p, (r) satisfies 

[d\ + 2'KdQ + k'^ + a^V{(t>ci)r]vn'^] ^/,fc(^) = . (AlO) 

1. Scalar perturbations 

Using Eq. ( IA8I ). we give the leading loop corrections to the correlation function of the stochastic variable ^ab- To compute 
the loop corrections to the scalar perturbations, first we have to find the correlation function of the density perturbation of the 
stochastic variable, which is given by the (0, 0, 0', 0') component of the noise kernel. It is obtained from the real part of i^ooo'O'- 
The part of order of {kH)'^ is 

F^'^Yo'iri, r2,k,p) ^ I d^x, I d^x,e-^f'-^^e-^P-^^F'^''^Y,,{x,, X2) 

= \ S^'\k +p) J d'q{{a,)-^[d,p,^'-^-q- {k - q)) + afVeU«'} 

x[ia2)-^[d,,'d,^'-'^-q-ik~q)) +4'Vci,2«'} G+(ri, r2)G+^_g|(ri, t^) 

(All) 

We put the momentum superscript on the partial derivative operator to represent that 9^ ' operates only to the Wightman function, 

G+(ri, T2). For example, d^^d^^''^ G+(ri, 'r2)G|^_^|(ri, T2) means G+(ti, '^2)9riG|^_^|(ri, T2). Taking into 
account that the correlation function of the stochastic variable is given by the noise kernel (12.21 ). we find Eq. (lAl II ) implies 



(<5p^fc(ri)5p,p(r2))W ^ j d'x, J d^x,e-^k-^^e-^P-^-^{e,{x,)e'A^,)Y'^ 



F(^)°„°;,(n,r2,fc,p)+F(^)°„°;,(n, r2, -fc, -p)' 



i<5(3)(fc+p) j d^q[{a,)-^(d,^d,^-^-q-{k-q))+afW,,,,K^] 

y.[{a2)-^(d,^d,^-''~q-{k-q)) + ^'V.z.aA^^jRe [ G+(ti, T2)G|^_^|(ri, T2) 



(A12) 
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On the third equality, we have transformed the momentum q to —q in the integral of f'^^^'^q '^qi{ti, T2, — fc, — p) * . This integral 
corresponds to the integral of the momentum of the internal line of the loop graph. 



2. Tensor perturbations 



Next we calculate the correlation function of the transverse traceless part of the anisotropic pressure for which we have to 
compute the loop corrections to the tensor perturbations. The correlation function of the pressure part of the stochastic variable 

is given by the bi-tensor i^'^'* in momentum space as 



Uri, r,,k,p) = I d'x.j d'x,e-^^-^^e-'P-^^F^'^) x,) 



(aiaa) " 5{k +p) I £q [{q'{kj - q^) + (F - q')}{q'' {k^, ~ q^,) + q„,'{k^' - q^' )} 



(trace part) ] G+(ri, T2)G+ _ (n, ra) 



(A13) 



The transverse traceless part of any tensor on spatial flat hypersurface can be extracted by means of the projection operator, 
Pj (fc), as follows. 



1 



= P\{k)P){k) - -P){k)P',ik) t\{k) 



P)(k)^6) 



k kj 



(A14) 



Operating the projection operator onto -F^*''* (n , r2, we extract its tensor part as follows: 

\F"'i '„,(n. n,k,p)]TT = 2(J^ +P) / G+(t,, T2)G+j._^|(n, n) 

X {-2,',, + - <?jf):) - + (i^) + 2^(.-,, + ,'*,)} 

(A15) 

Consequently, taking into account the definition of the noise kernel i2.2i . we can give the correlation function of the transverse 
traceless part of the anisotropic pressure pr^^*'' by 



(?'7r(*^(ri)e^(fc)p^J^(r2)e'',„,(p))W 



-{aia^r' 5(k+p) I d'q Re[G+^iT,, T2)G+^_^^{t^, 



{q-kf 



/s2 ) fc2 
{q-kf\ k^'k, 



q ■ k 



(A16) 



Especially when we contract the suffices (i, m') and (j, V), this correlation function is rewritten as 



(p7r^fc(Ti)e^(fc) pTT^p{T2)e\{p)) 



(4) 



\{a,a2)-^5{k+p) I d\ (9'-^ 



|2x2 



Re 



G+iT,, r2)G|^_^|(Ti, r2) 



(A17) 
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Appendix B: Loop integration 



To integrate over the inner momentum, q, it is convenient to consider the functions, f{k,S,J{) and g{k,5,'K) which are 
defined as 



f{k,5,:K) = 

g(fc,<5,J{) = 



qe[kiM] 



1 



1 



,3+5 



3+5 



\k-q 
{k-qf\^ 1 



(Bl) 
(B2) 



where 5 is an arbitrary small constant. Here we assume that ki is much smaller than fc, i.e., ki <^ k. In fact, fc,; is defined by the 
horizon scale on the initial time as ki = Jii. Since we are interested only in the modes whose scales are much smaller than the 
initial Hubble horizon scale, it is appropriate to assume ki <C k. 
To make the integration simple, we approximate \k — q\^(^+^) as 



for q < k 
for k < q 



1 



_L{l + 3tl + 0((,ffl')} 



\k-q\^+^ 
1 



Then f{k, 6, J£) is given by 



/(fc,<5,J£) = An 



1 1 



477 
fc3 



1 1 



l3 V/fc 



/ ^ 3 V5{ 



-2<5 



Similarly, g{k,S, J-C) is given by 



39Tr 

ff(fc,<5,J{) - ^ 
15 

327r 

Is" 



1 1 



{(fc - k^r^ - kf-'} + Y^s{^'-'' - (fc + k.r''} 



-2i5 3^_2i5 



(B3) 



(B4) 



(B5) 
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